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Abstract 

Agreement between the experimental value T{Ks — > 77) and the 
number predicted via a one-loop chiral perturbation theory calculation 
has been cited as a success for the latter. On the other hand the 
one-loop prediction for the closely related process ir^jj has 

been found to be a factor three below the experimental value. Using 
the inputs of unitarity and dispersion relations, we demonstrate the 
importance of higher order loop effects to both of these processes. 
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1 Introduction 



During the past decade, we have learned (at last) how to make rigorous con- 
tact between experimental low energy hadronic physics and QCD which is 
presumed to underlie such processes. This contact is provided by chiral per- 
turbation theory (ChPT) which exploits the (broken) chiral invariance 
of the light quark component of the QCD Lagrangian and provides a rep- 
resentation for interaction amplitudes as an expansion in energy-momentum 
divided by the chiral scale parameter 47rF^ ~ 1 GeV 0. A review of ChPT 
techniques will not be attempted here, but has been presented in a number 
of sources, wherein it is shown that, at least in the Goldstone boson sec- 
tor, such a chiral approach provides a remarkably successful and predictive 
representation of a wide variety of experimental processes 

The chiral technology begins by writing down an effective chiral La- 
grangian, the simplest (two-derivative) form of which is, in the Goldstone 
sector, 

4ff = -^TrD^UD^U^ + —TT2Bom{U + U^) + --- (1) 

where 

[/ = exp||l^^A,0,j (2) 

is a nonlinear function of the pseudoscalar fields, m = (m„, m^, ms)diag is the 
quark mass matrix, 

2B, = ^^ = ^^ = '-^ (3) 

TUu + ms mu + rrid rUu + nid + Arris 

is a phenomenological constant, is the covariant derivative, and F is the 
pion decay constant in the limit of chiral symmetry. Although these are only 
two of an infinite number of terms, already at this level there exists pre- 
dictive power-e.g., tree level evaluation of yields the familiar Weinberg 
predictions (at 0{p'^ ,m?)) for S-wave vr — tt scattering lengths ||^ which are 
approximately borne out experimentally. Loop diagrams required by unitar- 
ity produce terms of (9(p^,p^m^, rn^) and contain divergences. However, just 
as in QED such infinities can be absorbed into renormalizing phenomenolog- 
ical chiral couplings, and the most general "four- derivative" Lagrangian has 
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been given by Gasser and Leutwyler [|^ 




(4) 



Here the bare Li coefficients are themselves unphysical and are related to 
empirical quantities L^fi) measured at scale /i via 



where Fj are constants defined in ref. and e = 4 — d is the usual pa- 
rameter arising in dimensional regularization, with d being the number of 
dimensions. Gasser and Leutwyler have obtained empirical values for the 
phenomenological constants L^, . . . L^q. 

A wide range of electroweak and strong interactions of these Goldstone 
bosons have been successfully treated within this formalism, but there is at 
least one recurring problem — whenever the S-wave 1=0 tt — tt interaction is 
involved the simple one loop predictions have in general been found to be 
wanting This is perhaps not surprising, as the associated phase shift 
5o(s) passes through 90° somewhere in the vicinity of ^/s ~ 700 — 900 MeV, 
which has sometimes been associated with the existence of a broad resonance 
in this region. Such resonant behavior can certainly not be treated in any 
perturbative fashion and thus in this channel the chiral expansion must break 
down well before this energy is reached. 

There are a number of ways which have been used in order to avoid 
this difficulty. One is simply to confine predictions to a low enough energy 
that one- loop corrections should be sufficient. However, it is also possible 
to treat such effects nonperturbatively either by inclusion of S-wave 1=0 (cr 
or e) pole effects or by use of dispersion relations which effectively treat the 
TT — vr interaction to all loop orders — but at the price of introduction of model 
dependence This cost is usually considered worth paying, however, as 
by inclusion of such well-motivated model-dependence one can often reliably 
extend the usual region of validity of chiral predictions (i? <~ 400 — 500 
MeV) up to much higher values {E <~ 1 GeV). 

A particularly useful example of this procedure has recently been provided 
for the process 77 —>■ 7r*^7r^. Because only neutral particles are involved 
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(5) 
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there exists no tree level two-derivative or four-derivative contribution to this 
reaction, which guarantees that the one loop chiral perturbative prediction 
must itself be finite. This has been calculated as 

Here 



Fix,y) = 1 + ^ 

X 



1 -In^ ' 



X \p{s)-l 



with P{x) = J^^. (7) 



X 



However, the associated 77 — > nn cross section is given by 



-'^^ - - wi^-^r ' ^iF(v<)i^ (8) 

and is found to bear little resemblance to its recently measured experimental 
analog, as shown in Figure 1. 

The solution to this problem has recently been explored by a number of 
authors and has been found to be related to the inadequacy of the one-loop 
approach in the 1=0 S-wave vr — vr channel The solution is most clearly 
presented in terms of a dispersion relation approach. We assume, consistent 
with the chiral expansion, that when we are in the near-threshold region 
the only relevant higher order effects are in the helicity conserving S-wave 
amplitude, which we write as 

77-vr+7r-: f{s) = ^[2f,{s) + Ms)] 

77 -^V: r{s) = l[fois)-f2is)], (9) 

where I = 0, 2 refers to the isospin of the final vrvr state. For neutral pion 
production and working in the gauge wherein €2 ■ = £2 ■ ^1 = ci ■ ^2 = 
ei ■ fci = the transition amplitude is 

77^7rV: Amp = 2iehi ■ €2/^ (s). (10) 
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In the charged pion case the Bom and seagull contributions to this mul- 
tipole must also be included, so that the full amplitude becomes 



77 
Here 

where 



7r TT 



Amp = 2ze^ 



P+ ■ h 



ais) = 1 + fis) - fEUs 



-In 



I -Pis) 



.(11) 

(12) 
(13) 



is the Born approximation value for the helicity conserving S-wave multipole. 
In the threshold region the phase of fi{s) is required by unitarity to be 
equal to the corresponding vrvr phase shift 6j{s). When s > 16m^, inelastic 
reactions involving four pions are allowed. However, the inelasticity is small, 
being of order in the chiral expansion and will be neglected. 

The functions //(s) are analytic functions of s except for cuts along the 
positive and negative real axis. For positive s, the right hand cut extends 
from 4m^ < s < oo and is due to the s channel vrvr state. For negative s, the 
left hand cut is due to t, w-channel intermediate states such as 77r ^ tt — 77r 
or 77r ^ p ^ 77r, and extends from — cx) < s < 0. The single-channel 
final state unitarization problem has a simple solution in terms of the Omnes 
function 

e.p , 3 , ^4J(^] (X4) 



TT Jiml S' S' 



-the result must have the form 



fj{s) = gj{s)Dj\s) 



(15) 



where gi{s) is an analytic function with no cuts along the positive real 
axis. Morgan and Pennington consider a function pi{s) which has the 
same left hand singularity structure as fi{s), but which is real for s > 0. 
They then write a twice subtracted dispersion relation for the difference 
(//(s) -piis))Di{s), with the result lll| 



fjis)=Dj\s) 



pi{s)Di{s) + {ci + sdi) 



ds' pi{s')\mDi{s') 



7T Jiml S 



i2 



s — le 



(16) 



4 



where c/, dj are subtraction constants. Picking pi{s) to be given by its Born 
value and matching onto the known form of the low energy amplitude re- 
quired by chiral symmetry determines these constants unambiguously to be 



c/ = 0, di 



+ 



3847r2 



X 



/ 
/ 
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(17) 



One can now address the origin of the large corrections found in the 
77 — i> tt^tt" amplitude. Do they arise simply from the unitarization of the 
amplitude {i.e. Dj{s) ^ 1) or are new inputs needed in the amplitude? It 
turns out that the rescattering physics in Dj^{s) is most important, and that 
the main corrections are due to well-known ingredients. In our subsequent 
discussion, we will use a full phenomenological treatment but it is useful here 
to explore the case with a simple analytic form for Dj^{s). The condition 
ImDj{s) = —P{s)tf^{s) defines the [0,1] Pade approximation for the Omnes 
function |T^, i.e. 



with h{s) 



1 



1 - kis + t'f^{s){h{s) - h{0)) 
/?(s), fP{s) + l\ 



-In 



TT 



MO) = - 

71 



(18) 



and allows one an approximate but simple analytic representation for the 



77 



TTTT amplitude. The constant ko 



25ml 



is chosen to match the small 



s behavior of the experimental Dq ^{s) function, and ^2 = ~3o^ 
from a fit to / = 2 TTTT scattering. The resulting form for the 77 
amplitude is 



is chosen 



1 



;Fis,ml) 



X 



+ 



487r2F2 

i2s-ml)D,\s) + is-2ml)D,\s] 



3F2 



(19) 



which, when the Pade forms of Dj^^s) are used, provides a consistent an- 
alytic solution to the dispersion relation while also displaying the correct 
chiral properties to 0{s). In Figure 1, we plot the resulting cross section, in 
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comparison with the data and the lowest order result. It can be seen that 
the Omnes functions produce a substantial modification even near threshold. 
Of these, the most important is Dq^{s) which refiects the strong attractive 
TTvr scattering in the / = 0, J = channel 0. 

A much more satisfactory fit is found by use of an Omnes function Dq^{s) 
determined via the use of experimental values of the pion-pion phase shifts 
[0 as well as including contributions to the left hand cut due to Al,p,uj 
exchange diagrams, which leads to the very good fit given in Figure 1. Details 
of this calculation can be found in ref. [^. 

The lesson to be learned from this example is the importance of inclusion 
of 1=0 S-wave vr — vr rescattering corrections especially in processes which 
have no counterterm contributions and are generated from a simple one-loop 
chiral calculation. We shall see in the next sections two additional examples 
of this type, when we extend our formalism to cover the nonleptonic weak- 
radiative decays Ks 77 and — * 7r°77. 



2 K 



s 



77 



A good deal of work has been done extending the chiral formalism into the 
regime of nonleptonic weak processes. To two-derivative order the form of 
the effective SU(3) octet chiral Lagrangian is unique 



C 



AS=1 



FtG.TiiXaDMD^U^ 



(20) 



where is a constant whose value can be determined empirically. The 
corresponding four-derivative weak effective Lagrangian has also been written 



down and involves thirty-seven additional terms flj], whose explicit form will 
not be needed here, since as in the case of 77 — *■ h^tt^ there is no tree level 
contributions to the weak-radiative decay process Ks — > 77. Nevertheless 



there does exist a finite one-loop piece. Defining |]T5 



Amp (7^5 77) 



qi ■ q2 

the one-loop chiral prediction is found to be 



aF^B{ml) (21) 



m 



K. 



(22) 
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With the value 

G'i^ ^ 9.1 X IQ-^GeV-^ (23) 
determined from the tree level prediction for Kg — > tttt 

Amp\Ks vr+vr-) ^ 2F^A^^{mj,) = 2F^G^\ml - ml) (24) 

we find the one-loop chiral prediction to be in good agreement with the 
recently determined experimental number 

B.R.{Ks 77)chPT = 2.0 x 10^*^ vs. 
B.R.(i^s^77)e.p = (2.4±1.2) X 10-^ (25) 

An alternative derivation of this result is provided by the use of a disper- 
sion relation. We note that the function B(s) has a cut along the line element 



4m^ < s < oo and by unitarity has the imaginary part [T^ 

lmB{s) = 9{s - 4ml)P{s)Ao{s)f*{s), (26) 

where Aq{s) is the amplitude for the weak decay Ks tttt^I = 0) and fo{s) 
is the amplitude for S-wave radiative pion annihilation tttt^I = 0) — 77 
discussed in the previous section. Note that since we are using an SU(3) 
octet assumption for the weak transition, both di-pions are required to be in 
an isoscalar configuration. To lowest order in chiral perturbation theory we 
have 

A{s) ^ A^\s) = C^^is - ml) 

Ms) ^ W.) = ^^ln(^^-^j. (27) 

We now write a doubly subtracted dispersion relation for the function B{s) 
using as subtraction constants the requirements that the amplitude vanish at 
both s = and s = ml, as given in the lowest order chiral analysis, yielding 

^ s{s-ml) lmB{s')ds' 
^ ' TT 74m| s'{s' -ml){s' -s-ie)' ^ ^ 

If the lowest order chiral values — Eq. — are used to determine Im B{s') 
then the integration can be performed analytically to yield 

B{s) = -{ml~s)F{s,ml)G^\ (29) 

TT 
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which is precisely the (pion contribution to the) one-loop chiral result Eq. 
2|. 

At one level then the dispersive technique represents merely an alternative 
(and simpler!) way by which to perform the one- loop calculation. However, 
at a deeper level the dispersion relation provides a means to undertake a much 
more complete calculation of the radiative decay process by using for Im B{s') 
not just the lowest order chiral forms for these amplitudes but instead values 
which have more experimental validity. Of course, the subtraction constants 
must be fixed by some other means — we assume that chiral perturbation 
theory to one loop is accurate enough to describe the amplitude at very low 
energies. 

The Ks TTTT decay amplitude Ao{s) itself is an analytic function with 
a discontinuity along a cut 4m^ < s < oo given in terms of the S-wave 1=0 
TT — vr scattering phase shift 

ImAo(s) = e{s - Aml)e-<^'^ sin 5°(s)Ao(s). (30) 

The general solution of such an equation is given 

Aois) = P{s)D,\s), (31) 

where -P(s) is an arbitrary polynomial and Dq^{s) is the Omnes function 
discussed in section 1. In order to determine the polynomial P{s) we demand 
that the full Ks — > vrvr amplitude given in Eq. |31] match the simple chiral 
form Eq. in the absence of rescattering (i.e. when Dq^(s) = 1). We have 
then 

Aois) = Gs{s-ml)D^\.s) (32) 

while the corresponding solution for /o(s) was derived in ref. and has 
been outhned in section 1 of this paper.Q Substitution into Eq. |26| and nu- 
merical evaluation of the dispersive integral then provides a prediction for 
the Ks 77 amplitude which is much more complete and founded in exper- 
iment than is its simple one-loop chiral analog. The result of the numerical 
integration is shown in Table 1, where values are given for the radiative de- 
cay amplitude 5 (m|;) /branching ratio for four scenarios — a) simple one-loop 

-'^It is important to note here that when the solution Eq. ^ for Ao{s) is used, the 
value Gs = Gf^/|L'o"^(s = m^)| « 6.1 x 10"*'GeT^~^ must be employed in order that the 
proper normalization to the Ks tttt decay rate be preserved. 
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Input 


Im B [10-^] 


Re B [10-^] 


BR [10-6] 


ChPT 


0.77 


-0.44 


2.0 


Fade 


0.47 


-0.83 


2.3 


Gasser 


0.49 


-0.73 


2.0 


Gasser+V,A 


0.46 


-0.84 


2.3 



Table 1: Shown are the amphtudes and predicted branching ratios for the 
process Ks — >■ 77 using various theoretical inputs to the dispersion relation 
Eq. 28 as described in the text. 



chiral perturbation theory; b) use of the full unitarized values and with the 
Fade form for the Omnes function; c) use of the full unitarized values and 
with a numerical representation of the Omnes function based on experiment 
and given by Gasser [1^, 0; d) same as c) but with vector and axial- vector 
exchange contributions included also in /o(s), which provides the best fit to 
the 77 —>■ TTTT system. 

Examination of this table shows that the full unitarization procedure 
makes substantial changes in the predicted form of the decay amplitude with 
the importance of the real/imaginary pieces being interchanged for the ChFT 
and full unitarized calculations respectively. This is clearly seen in Figures 
2 where we show the very different shapes for Im,Re B{s) for these cases. 
Despite these differences the predicted branching ratio is remarkably robust, 
changing only slightly among the very different scenarios. We conclude then 
that the agreement between the experimental Ks — > 77 rate and that pre- 
dicted via one-loop ChFT should not be considered as a success for the latter, 
since we have produced nearly identical predictions for the branching ratio 
from a very different set of assumptions concerning the input parameters. 
Truong has reached a similar conclusion using a different parametrization of 
the Omnes function and an approximated 77 — > vrvr amplitude |T7 . 



Therefore, one must conclude that in order to distinguish between the var- 
ious decay mechanisms one must examine a process which offers the chance 
for a rather richer experimental confrontation — the related nonleptonic ra- 
diative decay Kl — > 71^77. 



9 



3 Kl^ 7r%7 

There has been considerable recent interest in the nonleptonic-radiative mode 
Kl 7r°77. This began when one-loop chiral perturbation theory was used 
to generate a supposedly reliable prediction 



Amp(K, ^ A7) = e^'e^ (^^^^^1^^) "^^"^ ^^^^ 



where 



C{s) = i-G^^ [F{s,ml){ml - s) - F{s,m],){ml + ml - s)] . (34) 

Since a three-body final state is involved, what emerges is a prediction 
for both the overall branching ratio 

B.R.{Kl 7r°77) = 0.68 x 10^^ (35) 

in addition to the shape of the decay spectrum, as shown in Figure 3. This 
distinctive shape is in marked contrast to that arising from a simple r], Vj pole 



model, which gives support at lower values of s^^ [^]. When experimental 
numbers were provided, the shape was found to be in good agreement with 
the ChPT prediction. However, the measured rate was nearly a factor of 
three larger than given in Eq. ^ 

^ / (1.7 ±0.3) X 10-6 NA31[20] 
B.R.{K, - vr 77) = I gg ± q g) ^ ^g-e fNAL|21|. ^^^^ 



Since this finding a number of authors have examined this problem. The 
inclusion of the A/ = 3/2 weak interaction results in a minor effect, as 
expected |22]. A dispersive analysis including unitarity corrections of Oijf'^ 
due to the vr+vr" intermediate state has also been presented recently [^. In 
this approach, the experimental results for the branching ratio as well as the 
spectrum in the invariant mass of the two photons can be reproduced if a 
somewhat sizeable contribution of vector meson exchange to the counterterm 
lagrangian of Oijf'^ is assumed. A similar result has been obtained in ref. 
p^ , however without taking into account unitarity corrections of 0{j)^\ 

Here we wish to examine the contribution of higher order diagrams to the 
decay process, by generalizing the dispersive approach which was applied in 
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the previous section to Ks 77. By definition only the amphtude propor- 
tional to the Lorentz structure in Eq. ^ enters the calculation. However, 
for this amplitude, unitarity corrections due to vrvr intermediate states will 
be treated to all orders in the chiral expansion. We begin by rederiving the 
one-loop ChPT result in this fashion. The amplitude C(s) possesses a cut 
along the real axis from 4m^ < s < 00 with a discontinuity determined via 



unitarity, which takes the form |^| 



ImC(.) = e{s - ^ml)\f5{s)A%Us)f\s) (37) 

where the lowest order ChPT amplitudes for 77 tt+tt" and tt+tt^tt^ 
are given respectively by 

f{s) = /^°^'^(3), 

Amp(ir^ ^ tt+tt-tt") = A^^,is) = G^^{s-ml). (38) 

As for Ks — > 77 we now write a twice subtracted dispersion relation for 
the function C(s), the subtraction constants being specified by the one- loop 
ChPT requirement that C{s) vanishes at both s = and s = m^: 

^ ^ s{s-ml) j-oo lmC{s')ds' 

vr J4ml s'{s' -ml){s' - s -ie)' ^ ' 

Using the lowest order chiral expression to determine Im C{s) the integration 
can again be performed analytically yielding 

C{s) = ^G<i^F{s,ml){s-ml) (40) 

which reproduces exactly the (pion loop contribution to the) one-loop ChPT 
result, i.e. the first term in Eq. 

To provide an improved estimate, we use the dispersion relation Eq. ^ 
but with a more complete representation of Im C{s) than just the lowest 
order ChPT expression (the subtraction constants are still taken from one- 
loop ChPT). In addition to the tt+tt^ intermediate state considered so far we 
include the discontinuity from the vr^vr^ intermediate state in the s-channel 
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to the unitarity relation 

(41) 



lmC{s) = 9{s - 4ml)l(3{s) 



where the superscript S indicates the S-wave component of these amphtudes. 
Other intermediate states open up at much higher thresholds and are sup- 
pressed in the twice subtracted dispersion relation. 

To evaluate the improved imaginary part of the function C(s) in Eq. ^ 
we shall employ dispersion relations to calculate the 77 — > vrvr scattering 
amplitudes and the — >• Svr decay amplitudes A+-o, ^ooo- The 

approach to has been reviewed in section 1. We shall use the explicit 

results given there. As for the ii" — Svr amplitude we shall use an approxi- 
mate solution to twice subtracted Khuri-Treiman equations Define the 
A/ = I 7^ — > Svr decay amplitude as 

km^{Kl ttVV") = 5'^*5"=^F(s„, Sb, s,) + permutations (42) 

where Si = [k — qi)^. Also note that F{sa, Sh, Sc) must be symmetric in its 
first two arguments according to Bose statistics. The decay modes relevant 
to our calculation are then 

Amp{Kl tt+tt-ttO) = A+_o = F(s„, s^, s^), 

Amp(K^ Sir^) = Aqoo = F(sa, Sf,, Sc) + permutations. (43) 

The discontinuity of F{sa, Sb, Sc) is provided by unitarity. Thus for K 3n 
the unitarity condition reads in general 

Im < 7rVV|7^||is:° >= 

algi (2;0^---(2;0^^'"^ ^(fc-P.-Pe-P/) 

X < TrVVlvrVvr-^ >*< TT'^n^n^nl\Kl > . (44) 

However, this equation is not amenable to an exact solution, so various ap- 
proximations are necessary. We begin by approximating the Svr scattering 



^This intermediate state does not contribute to Im C in the chiral one-loop analysis 
since vanishes at order . However the two neutral pions also have a / = component, 
therefore final state interactions in this channel are expected to give large corrections. 
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amplitude by a sum of two-pion, one spectator amplitudes — 

< ■K"'n''n''\7r'^n''n-^ >~ 5"-^ < ■K°'n''\'K'^'K'' > +permutations. (45) 

Also, since any two-pion reaction is at low energy we include only S-wave 
scattering terms 

+ (46) 

where 

A^'\s)^e<^'hmS'o{s). (47) 
The unitarity relation then reduces to the simpler form 

lmF{s,t,u) = l^[A('>(s)F{s,t,u) 

+ ^(AW*(s) - A^'>is))iF{u, s, t) + F{t, u, s)) 

+ ^A('>{t){F{s,t,u) + F{u,s,t)) 

+ ^A'^'>{u){F{s,t,u) + F{t,u,s))]. (48) 

The Khuri-Treiman equations arc obtained by the ansatz F(s, t, u) = U{s) + 
V(t) + V{u). We use a twice subtracted form with linear subtraction poly- 
nomials 

Uo{s) = au + bus, 

Vo{s) = 0. (49) 

The resulting system of equations has a simple approximate solution provided 
that we ignore the generally small 1 — 2 scattering term with respect to its 
much larger / = counterpart: 

V{s) = 0. (50) 



13 



[10-8] 


exp. 


Fit (Fade) 


Fit (Gasser) 




91.7±0.3 


89.6 + il9.6 


89.9 + il8.2 


/3i 


-25.7 ±0.3 


-21.8 - zl3.6 


-22.6 - zl2.2 



Table 2: Shown are values for the A/ = ^ component of the Kl 37r decay 
amplitude with final state interactions generated by the Pade and Gasser 
forms for D^^ respectively. 



where Si,S2 are the subtraction points, $2 is a correction which can be shown 
to be small, Dq^ is the Omnes function and [p^ 

2 r dVt 

U{s) = Uois) + -J —Uoitis,cose)). (51) 

We still need to specify the subtraction polynomial Uo{s). Since we want 
to input accurately the information provided by unitarity, we would like to 
use as much as possible the experimental information on — > decays. 
The K amplitude may be expanded around the center of the Dalitz 

plot as (neglecting A/ = 3/2 contributions) 

= «i - /3iY + (Ci + ^i)Y' + l{Ci- ^i)X' (52) 



with 



3 



Y = X = (53) 

mi mi 



Assuming the coefficients . . . are real, experiment fixes at the 

values given in column one of Table 2. As is well known, using the current 
algebra expression Uq{s) = Gg^{s — m^) as a subtraction polynomial gives 
values for ai and (3i which are too small. On the other hand, ChPT to 
one-loop can fit the K ^ "i-K data by adjusting the counterterm coupling 
constants of order p^. In the language of dispersion relations, this means 
that beyond leading order ChPT not only final state interactions contribute 
to coefficients ai, j3i but also the subtraction polynomial Uq is subject to 
corrections. This might also include contributions from higher resonance 
exchange. 
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In the same spirit we adopt the following phenomenological procedure. 
The K ^ 3n decay amplitude is calculated according to the approximate 
solution of Khuri-Treiman equations, Eq. The subtraction constants au, 
hjj of Uq are treated as free adjusted such that the experimentally found ai, 
(3i are reproduced. The values needed for the case of the Pade form of Dq^ 
are 

a^' = -2.3 ■ 10"^ = 8.6 ■ lO^^GeV"^ (54) 

to be compared with the current algebra expressions a^r^ = —Gsiri^ = —1.8 ■ 
10-7, 6gA = Gs = 9.1 ■ lO-'^GeV-^ [] This yields the resuhs given in the 
second and third columns of Table 2 for two parameterizations of the Omnes 
function, the Pade form, Eq. |1^, and the numerical parameterization given 
by Gasser respectively. Interestingly, /3i develops a rather large imaginary 
part as can be seen in the second and third column of Table 2. 

The advantage of this approach is twofold. First we use the experimen- 
tally available information on coefficients ai, (3i - the shortcomings of a too 
low K ^ 3tt amplitude from soft pion theorems are thus avoided. Secondly, 
it provides us with a representation of real and imaginary part of the decay 
amplitude outside the physical region. This is exactly what is needed in or- 



der to make Im C(s), Eq. 41, approximately real. If instead we would use 

the experimental expansion of q, Eq. as an extrapolation, Im C{s) 

would develop an unacceptable large imaginary part at rather low energies, 
i.e. for ^/s > 450MeV. Since our solution to the K 3tt amplitude is sub- 
ject to several approximations, we cannot hope that the imaginary part of Im 
C{s) cancels completely. However, it cancels to a large extent, i.e. it never 
reaches ten percent of the real part in the region from two pion threshold up 
to ^/s^ 600 MeV. 

We are now ready to calculate the function C{s) using the improved rep- 
resentation of its imaginary part. We used two parameterizations of the 
Omnes function Dq^, the Pade solution Eq. ^ and the numerical repre- 
sentation given by Gasser. In Figure 4 a) the improved imaginary part 
Im C{s) is compared to the lowest order approximation, and in Figure 4 
b) the corresponding real parts calculated from the dispersion integral are 
shown. In the dispersive approach, both imaginary and real part of C{s) 

■^In a similar approach, ref. p8| , Truong obtains a reasonable representation of the 
K ~> St: decay amplitude using current algebra constraints for subtraction constants and 
appending a 'p-pole' contribution. 
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input 


BR[10-^] 




exp 


1.7±0.3 




ChPT (pion-loop) 


0.59 


0.39 


Pade-charged 


0.93 


0.68 


Pade-charged+neutral 


1.12 


0.84 


Gasser-charged 


0.93 


0.68 


Gasser-charged+neutral 


1.08 


0.81 



Table 3: Shown are the calculated branching ratios for Kl 7r°77 under 
various input assumptions. The second column indicates the contribution to 
the branching ratio from the absorptive part of the amplitude. 



are rather enhanced, already just above the two pion threshold. Calculating 
BR{Kl — >■ 7r°77) with these inputs we find a net enhancement over the one- 
loop ChPT result by a factor ^ 1.9. The results are summarized in Table 3 
where the branching ratio is given for several scenarios. 

The observed enhancement is clearly seen to come from two effects:^ 

i) The use of a corrected K ^ ?)Ti amplitude in agreement with data 
gives an enhancement factor of ~ 1.6 in the branching ratio. This has 
been noted earlier in a very simple approach one could just 
scale Gs in Eq- il to reproduce the experimental — ^ 37r amplitude, 
leading to a similar enhancement factor for BR{Kl — > 7r°77). However, 
a consistent calculation clearly should explain all relevant processes, 
K 27T, K 37r, 77 — tttt. Kg — > 77 and — > 7r°77 by the same 
method. We have explained in detail above how this can be achieved. 

ii) The inclusion of the neutral two pion intermediate state in the uni- 
tarity relation brings an additional enhancement factor of ~ 1.2. As 
mentioned before, this contribution is missed in all approaches where 
the TTTT — >■ 77 vertex is treated only at the Born level. Although this 

^It should be noted that our result with inclusion of charged intermediate state only is 
consistent with that previously calculated by Ko and Rosner using a simple one-loop 
approximation. We do not understand the discrepancy with the dispersive calculation of 
Truong who also includes only the charged intermediate state and finds a branching 
ratio of 1.3 x 10"^ 
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effect is moderate, it goes in the right direction. Moreover, it is critical 
to include the tt^tt^ intermediate state in order to properly implement 
the constraints of unitarity. 

Finally we note that using a twice subtracted dispersion relation, the 
contributions from the high energy region to the dispersive integral are very 
much suppressed. Cutting the integral at ^/s — 600MeV instead of 1 GeV 
changes the branching ratio by less than one percent. Also, the calculated 
spectrum in the invariant mass of the two photons is plotted in Figure 3. 
It deviates insignificantly from the spectrum obtained in one loop Chiral 
perturbation theory. If only the charged intermediate state is included, the 
maximum of is shifted toward higher (f values, contrary to the exper- 
imental trend. Inclusion of the neutral intermediate state in the unitarity 
relation restores the maximum to its original location. 

4 Conclusions 

The process Kl — > 7r°77 has traditionally been a difficult one to understand 
within the context of chiral perturbation theory. Indeed, from other suc- 
cesses one might have expected the one-loop chiral prediction to be accurate 
to (9(20%) or so, whereas in fact the predicted branching ratio is nearly a 
factor of three too small. Previous work in this area has attempted to ex- 
plain this discrepancy in terms of vector meson pole contributions and/or 
in terms of higher loop final state interaction effects. We have here noted, 
however, that previous dispersion-based final state interaction calculations 
have included only the intermediate tt+tt^ intermediate state, omitting its 
potentially important tt^tt^ analog. Above we have given a mutually consis- 
tent analysis of both Kl — > Stt and Kl — > 7r°77 processes using dispersion 
relations and have shown that by including this previously neglected tt^tt^ 
intermediate state piece the branching ratio for Kl 71^77 is significantly 
enhanced, although it remains too low to fully explain the data. This should 
by no means be considered to be complete analysis — indeed many effects 
such as 1=1,2 TTTT scattering effects as well as higher order contributions to 
the Kl — > Stt process have been neglected. Nevertheless, we believe that this 
calculation opens up interesting questions for future study about the impor- 
tance of effects beyond the one loop approximation in chiral perturbation 
theory. 
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Figure captions 

Fig. 1: Indicated are experimental data points for 77 — > nn compared to the 
one-loop cliiral perturbative prediction (dashed) and dispersive calcu- 
lations using Pade (dotted) and Gasser (solid) Dq^{s) functions. 

Fig. 2: Shown are the results of one-loop chiral perturbation theory (dashed) 
and dispersive analysis scenario d) (solid) for Im B{s) (a) and Re B{s) 
(b) respectively. 

Fig. 3: Normalized spectra l/TdT/dz in the invariant mass of the two photons 
{z = q^/Mj^u). Plotted are chiral perturbation theory (dashed) and 
dispersive analysis using Gasser Dq^{s) function (solid). 

Fig. 4: Im C{s) (a) and Re C{s) (b) in the physical decay region. Shown are the 
results of one-loop chiral perturbation theory (dashed) and dispersive 
analysis using Gasser Dq^[s) function (solid). Also shown is the result 
of the dispersive analysis including the charged intermediate state only 
(dot-dashed) . 
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